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ABSTRACT 


The present work is based on simulation ot boiling 
by air injection into water medium to obtain the 

frequency response of bubbles for air flow rata range 
128.9 cc/sec to 644.5 cc/sec and water flow rate range 
321,88 cc/sec to 431,85 cc/sec. 

During the bubble collapse the acoustic signals are 
picked up by a piezo electric transducer. The signals are 
amplified/ filtered and then it is auto— correlated. To get 
frequency response, the auto-correlation is fevd into 
the Fast Fourier Transform Analyser. The power spectra 
obtained for the system contain information about the bubble- 
size, shifting frequency as a function of velocity of 
the liquid. 

The dependence of predominant bubble frequency is found 
to be directly proportional to the 4-th power of flow 
velocity. 



CHAPTER I 


INTRODUCTION 


The term noise when applied to reactor sti-'dy/ it 
means random phenomena like the fluctuations in neutron 
population. Analysis of noise sources in a power reactor 
can give vital intormation about the state of the reactor 
without disturbing the normal operation. According to 
Fourier's theorem, a function of time is capable of being 
cnarac tarized by tne superposition of sine waves of various 
frequencies having appropriate amplitudes and paases. The 
amplitudes that are present at various frequencies therefore 
become an alternative method of characterising the noise. 

In general a random variable x, which is a function 
of time can be represented as 

CO 

x(t) = / X(f) exp(ift)df 

-OD 

where X(f)df is the complex amplitude of the oscillation 

in the region between f and f+df. The behaviour of 

2 

magnitude/ lX(f) I , with frequency is one of the most useful 
methods of cnaracterising the noise. 

Some of th i applications of noise analysis in 
reactors are power monitoring,, coolant flow— rate and void 
fraction measurements, det-action of nucliata boiling and 
structural vibrations. The same noise analysis techniques are 
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used in industries. Where a large capital investment is 
involved in a plant, it is expensive to stop production 
in order to do testing. Hence noise analysis is a powerful 
tool for monitoring a system during its on-line process. 

The objectives of the present work are: 

(1) To demonstrate that noise techniques can be 
used to obtain information about nucleate boiling, 

( 2) To show that noise techniques can be used to 
dv:;tarmine the frequency response (spectral density) of 
bubbles in nucleate boiling for various flowrates of 
coolant. 

The present work is bashed on siiaulation of boiling 
by air injection into water medium and to obtain the 
frequency response of bubbles. The power spectra obtained 
for the system contain information about the bubble— size, 
shifting frequency as a function of velocity of the 
coolant. 

One of the principal advantages of the correlation 
techniques w-rxich is used in this experiment is to extract 
signals buri x3 in noise. 

The experimental equipment and procedure are descrdbad 


in Chapter ( 7 ) 


CHAPTER II 


EXPERISHCES ICEIH NOISE AICALYSIS 


2.1 Exp^.rience in Lo«se Parts f-'Ionitoring and Operating 

Nuclear Power Plants 

Some of the earlier investigations are given below; 

The monitoring of the french "PHENIX" LI4F3R [1 ] 
dynamic behaviour is assumed by periodical recordings of 
signal fluctuations nf neutronic noises and vibration 
measureiTients. The frequency analysis of these recordings 
made by the power spectral density and the coherence function 
analysis allowed the interpretation and monitoring of most 
noise sources caused by pump functioning and structure 
movements under sodium flow. The results that are obtained 
are as f ol lows : 

(1) A movement of the concrete slab of the primary 
circuit at the fundamental frequency corresponding to the 
pumps rotary speed . 

( 2) A low level vibration on a control rod drive 
corresponding to an incipient defect. 

(3) Emergence of vibrations on the core cover 
during a working cycle. 

(4) A resonance corresponding to the first vibra- 
tional mode of the fuel assemblies under the cooling sodium 


f low 
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The evidence of the vibratory movement has been 
obtained by the systematic calculation cf the coherence 
function between the available neutronic noise measure- 
ments and the vibratory transducers. The experiences 
gained on PHENIX vibration measurements are included in 
the SUPERPHEin;x monitoring system. 

Fry and Robinson [2] observed a sudden peak in 
their library of neutronic noise soectra. The peak was 
attributed to the control rod bearing failure. The peak 
disappeared when the defect was removed. Hence the above 
example illustrates now noise spectra were used to find 
the mechanical failure of a system. 

2,2 Boiling Detection in Fast Reactors by Hoise Analysis . 

Studies Performed in France 

The onset of accidental conditions particularly 
the boiling of sodium is generally related to background 
noise modifications in an operating reactor. Tide detection 
of in-core boiling may be achieved by three different 
analysis methods. They are: (1) neutron noise, C 2) acoustic 
noise and (3) temperature fluctuations at the subasssribly 
outlets. On these three methods experiments ware conducted 
in RAPSOBIE reactor in France, These results are as follows [3], 

Thermohydraulic analysis of boiling in CFNa bundle 
have demonstrated that two main regimes may appears 
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- tint spot boiling 

- single buoble ragima. 

Tha cransition between these two regimes is fairly 
sharp. Tho behaviour of tha boiling process is rather 
insensitive to thermo hydraulic parameters (flow/ heat flu3<^ 
geometry) , 

It appears that 

— Thermal noise detection is especially effective 
in hot spot boiling conditions* 

— Acoustic noise is directly sensitive to hot spot 
boiling and also, modulat-ad by tha single bubble regime, 

-• Void effect maasures either by pressure or flow 
fluctuations or neutron power fluctuations, is only detected 
in the presenc.' of single bubble. 

The above mentioned three methods ara complementary 
and correlation between than leads bo a reliable boiling 
detection system, 

2 ,3 Literature Survey 

Strasberg [4 ] has given that gas bubbles when 
entrained in water or other liquid can generate high sound 
pressures in the liquid. Significant sound pressures are 
associated only with \7olume pulsations of tho bubble whereas 
oscillations in the shape of the bubble do not result in 
appreciable sound. He has done calculations of sound 
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pressures resulting from excitation of volume pulsations 
by the following mechanisms: 

(1) bubble coalescence, or division? 

( ii) The motion of a free stream of liquid containing 
entrained buboles past an obstacle? and by the flow of 
liquid vjith entrained bubbles through a pipe past a cons- 
triction. The calculation of the sound pressure generated 
by bubble formation has been verified by measurements with 
bubbles formed at a nozzle, 

Claytor [5] has injected steam into sodium to measure 
their acoustic cnaracteristics. Steam was injected through 
small swaged leaks in Type 304 stainless steel tubing. The 
swaged portion of the leak was typically 2.5 cm long. During 
a leak, acoustic noise was measured simultaneously with lithium 
niobate microphone end accelerom.^ters. He repeated with 
inert gas injection into sodium. And concluded that 
frequency peaks upto 20 KHz are due to gas bubble resonance. 

^it frequencies between 20 and lOO KHz, combustion noise is 
probably the main source of noise for steam injections. Above 
120 KHz the mechanism of sound generation by steam and 
inert gas injections are not fully understood. It is 
concluded that noise generation is dependent of sodium 
cover gas pressure, 

Claytor [ 6 ] also conducted an experiment at ANL for 
various sodium flow at various steara injection rates. He 
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concluded that background noisa at 2 KHz peak is due to 
argon hubbies at the orifice of the tube which was to keap 
the injector from plugging. An increase of 8,5 dB over 
the background level is observed at 2 KHz when steam is 
injected. The noise produced at 2 KHz was possibly 
attribute to acoustic monopole radiation from hydrogen 
bubbles oscillating in the volume mode aftsr being formed by 
the steam-sodium reactions. Finally, it was snown tnat 
bubble-ganerat^ad sound due to peaks ara prevalent at low 
frequencies ( f < 10 KHz) as is the ambient noise which is 
present in a LMFBR— SG du a to flow^ 

Ying [7] developed an acoustic detection system to 
pr av'ant serious damage to the tubes that would result from 
sodium water reactions. The leak can be identified from 
the acoustic spectrum of noise generated during sodium- 
water reactions which produce hydrogen gas with sodium hydroxide 
or sooium monoxide. The significant audible sound is the 
sound radiated from tha damped oscillations of hydrogen 
bubbl as in liquid and the frying noise generated at the 
interface of sodium and water have peaks in tha acoustic 
spectra at frequencies below 1,5 KHz and near 2 KHz 
respectively '"'^u a to exothermic heat of the chemical reactions, 

DB [8] measured the dis tribution of time intervals 
between incipient bubbles in a venturi and concluded from 
the observations that cavitation noise observed is neither 
the result of an impulsive random pulse nor a purely cyclic 
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process. Such noise results from bubble clusters incepted 
periodically in a train ( 3 to 5 appearances) and trains 
appearing at random. This is possibly due to the back 
action of bubble nucleation on flow turbulence. The above 
observations provide further understanding of the fundamental 
process and preliminary data that indicate possibilities 
for detection schemes in monitoring systems for incipient 
two phase flow, 

Gavrilov [9] showed the results of measurements of 
the sound attenuation induced by gas bubbles can be used to 
find the size distribution of the bubbles in a long-standing 


water and tap water 



CHAPTER III 


SIGNAL PROCESS li^G 


Th<= most simple descriptor used for getting informa- 
tion from fluctuating data is the 'mean' value that expresses 
the steady f low of information in the data or signal. 

To analyse the signals coming from more complex 
systems and the limitations of the simple mean and/or 
variance parameters led to the evolvanent of better statis- 
tical descriprors such as correlation functions, power 
spectral density, probability density distribution etc. For 
these calculations we require analog ani digital correlators, 
spectrum analysers based on Fast lourier Transform Techniques, 

The use of fast Eourier Transform Techniques developed 
in mid 1960's has bean really most revolutionising and many 
of the tasks that were possible only off-line before their 
advent/ can be periormed on-line today in the plants. 
Development of digital filters, fast analog to digital 
convertor and the use of multiplexer have been useful in 
studying the several signals simultaneously along with their 
cross correlation. 

The availability of analog tape recorders nave made 
storage of data easy and hence the scope of extensive 
research on available signals. 



CHAPTER IV 


METHODS CF MEASURB.MEI^IT 


There are two mai-hods of noise recording and 
analysis; 

(1) Analog (or continous)/ and 

(2) Digital .or discrete). 


Continuous recording may be accomplished by standard 
chart recorders or by magnetic tape recording. Having 
once recorded an exceriitent on magnetic tape, it can be 
return and reanalysed at any time and as or ten as is desired. 
Data so recorded is easily stored and is compatible with 
both analog and digital methods of aiialysis. Also by using 
different recording and playback speeds/ one can achieve 
various time transformations in the analysis. Although 
continuous signals can be recorded such that the degree of 
magnetization is proportional to the signal (direct amplitude 
recording) / it is usually preferable to record the frequency 
modulation ot a carrier at constant amplitude. Auto— correlation 
Functions can be computed from the playbacK-head signals 
rom identical channels by varying the length of tape between 
le two heads to achieve various delays. 


It is vary common in reactor installations to have 
? reactor power and other plant variables recorded on X-Y 
’tters. The charts give the basis for off-line digitizing i.e. 


r reactor operation. It provides a permanent record of 
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the experiment. Experimentalists may v.’ish to ‘sae’ what is 
nappening during the experiment. 

When digital analysis of a continuous signal is 
desired/ an ana log- to-di gi tal conversion muse be performed. 

The schematic diagram of signal processing is given 
in Chapter ( 7 ). This is an on-line electronic analysis 
i.e. during plant operation. 



c:~:\PTSR V 


IvIATHgMATICAL Ai'JALYdIa OB 


■ 5 . 1 Theory 

Variance 

The brief treatment of random signals (noise analysis) 
vj'hich is given in this section follows essentially that of 
Thie [10] and Rice [H]. For a random process^- the n-th 
moment having proPability density function p(x), is given 
by 


x"" P(x) dx 


/ p( x) jx 
-cp 


( 1 ) 


or, for a discrete distribution. 


S x" D(xp 
k '’<^k> 


( 2) 


Equation (2) suggests that x^ may be determined 
experimentally from a sequence of a large number, of data 
values, Xj_,/ since these will tend to distribute themselves 
according to p(X|^): 


n 


X 


i'^ 


N 

2 

k;=l 



/ 


or, for continuous data. 
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X 


1 

T 


T/2 

/ 

-T/2 


dt. 


whora T is the duration of tha experim-ant for wnich x( t) 
is available. It is avid ant tnat the mcmants are average 
values of various powers of tha random variaole. 

The first rncmant, x^ is tha average value, for most 
noise applications it is convenient to choose the co-ordinates 
such that X will be zero. But to avoid loss of generality, 
it need not be zero here. That being the casi, it is 
more convenient to avaluata moments with respect to x. 


^ (x-x) p(x)dx 

QO 

I p(x) dx 

CO 

For n = 2, the result is called the "variance". The square 
root of the variance is termed tha standard deviation h : 


, n 
(x-x) 


2 

a 


GO __ 2 

/ (x-x) p(x)dx 

-oo 

CO 

p(x) dx 


1 

T 


T/2 __ 2 

f (x-x) dt 

-T/2 


For digital data, 

N 

2 2 
a lc=l 


S ( b( Xj^) 


N 

2 P(x.> 

k=l 


1 

N 


N 


k=l 


(x,^ - x) 
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NoW/ 


1 

N 


N 

S 

k=l 



1 

N 


/ 2 - - 2 , 

2 (X.-2X X + X ) 

k=l ^ ^ 


N XjJ 

2 — 

^ N 

k=l 


2x 


N X 

2 IT + 

k=l 


N 

2 

k=l 


-2 


X 


N 



0-2 
2x 4 - 


-2 

X 


2 


c 



Tner^fore standard :i aviation. 


O' = 




1/2 


The standard deviation is also referred to as the 
"rms" (root mean square) value. The rms value is very 
commonly used as a quantitative measure of the amount of 
noise, since it is easy to compute, it has particular signi- 
ficance in those distributions that are gaussian, and it is 
easy to measure exnerimentally. 


5.2 Correlation Functions and Prv;er Sp>actral Density Functions 
We consider a physical process which gives rise to a 
time-varying signal x(t)* The signal inay oe simple or 
complex periodic or have the character of noise, that is 
random- varying, Ir tnis signal is delayed, thereby obtaining 
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the signal y.(t-T)/ whicn is identical to xi. c ) buc just 
delayed in time, and then the product o± x(t;«x(t-'t) is 
averaged over a sufficiently long time, the auto— correlation 
function, 0 ('c), for the signal x( t) will be determined. 

Mathematically, the auto-correlation function is defined 
as. 




Lim 
T -* oo 


2T 


T 

S 


-T 


x( t) xCt-T) dt 


( 1 ) 


Note also that 0 (f) is an even function, that is 

= 0^^('-'^), anr therefore eqn.(l) may be written as 


= Lim 


rn 

X 


CO 


T 

/ x( t) x( t+'t) dt 
-T 


( 2 ) 


For digital signals having N discrete data points, 
eqn«( 2) can be re-written as replacing integral by summation 




Lim 


w ^ oo 


N 

^ k=l 


U) 


The following are the properties cf this function. 


0( t) ; 


1 . 

This states that the Auto-correlation function is an 


even function and is symmetrical about T = O 
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2. = mean square value 

This states tnat ac t = O the total signil pov;er 
(AC and DC) is represented. 


3. 0^„(oo) = (average value) 

This states that the value for large values of t is 
approaching the DC power of thd signal. 


4. 


1 

27f 


CO 


, , / i 2'iif . ^ 

J § if } e or 


-CD 


This states tnat the Auto-correlation function and 
the power spectral density §(f) forin Fourier Transform 
pairs. 

5. In Auto-corr elat-cn calculations, pnase information is 
lost. 


The normalized auto- correlation function is defined 
w.i.t. to deviation from the mean and is given by 


P^jj 5 ^(t) = Lim 

T -» CO 


2 cr^ T 


T _ 

/ p( t)-X][ x( t+T) - X ] dt 

-T 

(4) 


lor digital techniques. 


= Lim 

N -» CO 



N _ 

2 [X ( tj^)-x ].^x( tj^+f) -X ] 
k=l 


where x is the mean value of the tims varying signal over 
one period of a periodic signal or over the discrete number of 



17 


2 

digitized values of a random signaly and a is tine variance 


ot the signal x( t) . In practice, with N data poiiits spaced 


h t apart, is calculated for t upto t = MAt <i''I At 

ixx m 

where M is the maximum number of correlation points. 

“ is the total time delay span which is equal to M t. 


“ra 


0 (^) 

XX 


wtaere x( t) 
and 

2 

a 

XX 


1 


i<I- T/At 


N-V'-^t k:=l 




( 5 ) 


XX 


1 


1 

,J 


M 

S t) 

K=1 


N 

k=l 


( X, 


2 

x) 


The auto— correlation function and the normalized auto- 
correlation are related by 


0 (t ) = 0 (t) + Cx)^ 

^XX XX XX 


( 6 ) 


The correlation functions are useful in describing a 
system's response in the time-domain. 

Now, the frequency response of a system cun be 
obtained by talcing the Fourier transform of the correlation 
functions. Then useful reciprocal relations are known 
as Weiner's theorem. Namely, the power density spectrum, 
^(f), of a signal is the cosine Fourier transform of the 
auto-correla tion function, 0 (t), 



Tnis may be expressed as 


PC') 


CD 

/ $ ( f ) exp( iwt) df/ w = 2Jv£ 

-CD 


CO 

= 2 / 
o 

where^ 

^(f) = Lim 

T -♦ CD 

§ i±) = Lim 

-* CD 

m 

= Lim 

It is evident that 

= x^ = 


f(f) coswT:df 


^ lx(f) i ^ 

■"m 

/ 0{x) expC - iwT ) d x 



2 j 0{x) cosw X dx. 
o 


CD 

/ § { f ) df 

-oo 


( 7 ) 


( 8 ) 


2 — 2 

is the total power^ v;hich by virtue of a = 0(o) - (x; is 

2 - 2 
made up ci an a-c. Component a and a d.c. component (x) . 

The usefulness of the auto-correlation function therefor 

extends beyond its presentation of information in the time 

domain. By Fourier analysis the correlation function can give 

information in the frequency domain via the power spectrum. 
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^ ^ Cross-Corr alation i?‘unction and its Properties 

tv70 random functions x and y are monitored, the 
cros s-corr 3lat.ion is defined by the folioi;ing equations: 

1 ^ 

= him — / x( t) y(t+T^)dt, 

^y IT 

T ^ oo 

1 

^yx(^) 2 T Y<. t) x(t+T:)dt. 

T ^OD 


This tormula is a point-by-point multiplication of a 
waveform by the shifted second waveform, followed by an 
integration or summations over all time- 

The following are properties of this function; 


1.0’ C -T ) = 0 ( t) . 

^xy yx 

The cross-correlation displays symmetry about the 
ordinate when x and y are interchanged. 


2. I < . 0yy(O) (i) 


The equations define the bounding relationships for 
the magnitude of the cross-correlation function. The first 
states that the square of its magnitude is never greater 
than the product of the power contained in tne two signals* 
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The second s s that its magnitude is never greater tnan 
the average of the power contained in the two signals. 


^ ( f ) 

^ xy 




ph -5 eo'ove -.guations states tnat a Fourier Transform 
pair relatlc'iisiilp exists between the cross-correlation 
function ■■-n rhe cross power spectral density function. 
Unlike the Autc-correlation^ the cross-correlation contains 
phase informaticii besides amplitude. For digital techniques^ 
the cross-correlation function is 


$ 2 ^ 




M 

S =<( t;,) -!'< 


where M is the number of points 



CHAPTER VI 


ELBCTxRONxC IHSTRUMa.\^TS DESCRIPTION 


Piezo-elactric Transducer 

Transducer is a device which converts one form of 
energy into anotnar form energy. 

For example Icudspealcer is a transducer which 
converts electrical energy into acoustic ?nergy. Micropnone 
is a transducer whicn converts sound energy into electrical 
energy. 

In this axoerimenc a piezo— el ectric transducer which 
has a flat (+ 3 dB) frequency response up to 2 i'4HZ/ made by 
ECIL is used cc convert acoustic signal, caused by change in 
pressure when bubble coliapsas/ ‘into electrical signals, A 
Microphone which has a flat + 3 dB freqeincy rasponse 
upto 20 Khz, made by PHILIPS is also used to verity the 
signals that are obtained by ECIL transducer. 

6,2 Charge Amplifier (Model 27 35 ) 

The instrument was designed and manufactured by 
SNDEVCo, Division of Pec ton, DicKinson and Company, U.S.A. 

This charge amplifier is an all solid state instrument 
designed pacifically to condition the signals from a 
transducer. It will airplify the signals which are generated 
by transducer. 
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6,3 Tunable Band-Pass filter 

Band-pass filter means it will pass only certain 
frequencies and reject otner frequencies. 

In this experiment KOMBINATIONS FILT3R of type 01006 
and type 01007 ©f VEBRP'T MSSSELEKTRONIK > OTTOSCKN < DRESDEN^ 
West Germany is used as a high pass filter. It nas tunable 
cut-off frequency arrangements. In the present experiment 
a frequency of 125 Hz is used for nign pass filter. 

Correlation and Probability Analyser (Model SAI-48) 

The instrument vvas designed and manufactured by 
Signal Aria lysis Operation, Test Instruments Division, Honeywell, 
Denver, Colorado, U.S.A. 

The SAI-48 provides on line real tirre computation of 
Auto and cross correlation functions with incremental lag 
or time delays ranging from 0,05 [is to 2 seconds resulting 
in maximum time delays from 20|is ro 800 s^econds. An Auto 
or cross-correlation function is determined sequentially at 
400 incranental lag points so that a complete correlation 
function is displayed at one time. 

In addition it has two other primary operating iruades; 

1, Probability (Density and Distribution) 

2. Enhancement (or Signal Recovery). 

In all modes the 400 analysis points are computed, 
and may be displayed on external equipmenc such as an 
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osG il loscope^ or an X—Y recorder, 

6,4,1 Correlation 

Correlation analysis provided a quantitative measure 
of the deqree of similarity between waveforms as tney are 
-oeinq shifted relative to one another in time. If the signal 
is being compared with itself^ the resulting waveform is an 
Auto-correlation, Any two different waveforms may be compared 
via cross-correlation. 

The detection of signals in noise/ the determination 
of dynamic system errors/ automatic adjustment and control 
of processing plants,- localization of interference sources/ 
directional recBption of signals, determination of speech signals 
and evaluation of ballisto cardiograms are a few possible 
applications. Model SAI-48 is used to find auto correlation 
in this experiment, 

6.5 Fourier Transfoicm Analyser (Model SAI-470) 

The instrument was designed and manuf acrured by 
signal Analysis Operation, Test Instruments 'Jivision, Honeywell 
Inq,, Denver, Colorado, U.S.A. 

6,5.1. General Functional Description 

’ ■ "'\T r'l'"" ^ T'-' ■' '"'I |■|■|IIrll'T,-^■-rr^|^^ r nT 

The function of tne Model SAI— 470 is to compute the 
amplitudes of, and tne frequencies present, in an input function 
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presented to it via an SAI"48, 

The frequencies computed can be thomgnt of as frequency 
vectors Ci.a., complex frequencies) navinq a real and 
imaginary value/ and therefore/ a phase relationsnip. 

Most specifically, if represents the complex 

frequency magnitudes of some frequency then the SAI-47 lj 

solves the following basic equation: 

% kn 

M -i — - — 

^(£ )=S S2f (t) (e "'^ )/ (1) 

k=0 

where k = Number of samples representing the input tLme 
function and k = 0, 1, 2,3, . . . .l'I„ 

= magnitude of the k-th input of the input time 
varying function, i.e. correlated function 
(which is in time domain) 


“ i Txkn 
e N 


= cos 


nkn 
“n” 


i sin 


71 kn 

’TT" 


N = 1000 ^maximum number of frequency intervals) 
and n = 1,2,3,„„» 1000 . 

Furthermore eqn,(l) can be expanded to include its 
real, and its imaginary componsnts where 

^ ISOkn 

= kSo cosC—^) 
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and 3 


n 


M 

S 
k=0 


-i 2 sin( 180 


(3) 


From Gqns.( 2) and (3) the absolute fraquancy P 


n 


magnitude can be computed: 


n 


= [<V' (V"] 


211/2 


(4) 


and the angle can be computed from 


B 


tan / 


n 


n 


n 


(5) 


or 


sin 9^ 


n 


( 6 ) 


To signiticnatly smooth the power spectral density, 
the input function ) is multiplied by a Hamming weighting 

function. This weighting function is defined as 

m = 0,54 + 0.46 cos(- ~°- ) (7) 

where M = total number of data samples (M = 400) 
k = kth sample. 

Also, 1® multiplied by a scale factor of 1,28 to 

account for the fact that a maximum of 4oO summations occur, 

9 

The number 400 is closest to the binary number 512 ( 2 ) and 
therefore 512/400 = 1.28. This scale fcictor is actually taken 
into account in the weighting function as follows: 
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WF* = 1.28 [ 0.54 + 0.46 cos( --^ - ~ ) ] (S) 

or 

Wf’ = 0,6912 + 0.5888 cos 

I 

wuere WF is scale factored weightiiig function. 

Therefore, 

= 0^1.'^) [0.6912 + 0.58tib cosC ^- - ^- — ) ] (9) 

I 

0 (T) is the siTiootnened auto— corre lation function. 

The sinootnened auto-correlation function is used in Cl)- 

to find power sjiectral density. 
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IX?BRIMI,^ITAL BQUIPM£NT AlJD PR0C5DURS 


The experimental arrangement is explaine-i in the 
following paragraphs, 

A glass cylindrical chamber of size 4" diameter ant 
8" length is provided with an inlet and outlet. At the 
centre of the chamber a 2 m„m, hole is provided for air 
injection. In this chamber air is injected by means of 
a compressor at various air flow-rates. Now cnamber is 
filled with water. When air is allowed to pass through 
the chamber air bubbles are produced. When air bubbles 
collapse/ the acoustic signals are picked up by means 
a piezo electric transducer and converted into electric 
signals. These signals are amplified by the charge 
amplifier. 

The signals froa the charge amplifier for vario'^s 
air flow- rates are displayed on oscilloscope. These signc-.i 
are allowed to pass through a bind-pass filter and the 
lower cut-off frequency was set at 125 Hz, Thus the 50 ?. 
a,c, components ar-e eliminated. Now the signals that ar^.- 
coming from high pass filter due to bubble collapse are 
fed into signal correlatQT to compute auto correlation 
with delay time *5 m,s„ Now the correlated function is in 
time domain. The correlated function is converted intc 



frequency uomain by feeding tne output of the correlator 
intc a Fast Fourier Transform Analyser. 

The output of the i.F.T., magnitude versus frequency 
is displayed on the oscilloscope. From the scope, the 
frequency of thi acoustic signal due to bubble collapse is 
found. 

The experiment is repeated for various water f low- 
rat 3S (i. -3, various pressure inside the chamber) and air— 
flov^-rates. The pressure inside the chamber is measur-ed 
by means of a J-tube manometer using mercury’. 

The sch 3Tiatic diagram of the 
ment is shown in Figures (1) and (2) 


experimental arrange 
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CHA.PTER VIII 


EXPBRD4SNTAL OBSERVATIONS AiTO TA3UIATI0NS 


Air is injected by means of a compressor into the 
chamber at various air- flow- rates in stationary water. 

The power spectral cfensity at various air- flow-rates 
in stationary water are tabulated in Taole 1. 

TABLE 1 ; Frequency Spectrum of Air Injection into 
Stationary Water 


Air Flow Rate 
Gc/sec 

— 1— — - 

! with Bubbles 

! Frequency Mag, 

1 (Hertz) mV 

1 

~i 

J T/£Lthout Bubbles 

I Frequency Magnitude 

I ( Hertz ) mV 


153 

20 

153 

20 


203 

25 

203 

25 

128,9 cc/sec 

340 

30 




406 

40 




556 

20 





153 

20 

153 30 



202 

25 

20 2 35 


257.8 cc/sec 

320 

50 




410 

30 




5 53 

40 




198 

30 

201 20 



335 

20 



386,7 cc/sec 

408 

35 




543 

40 




201 

20 

2u2 20 


515,6 cc/sec 

331 

30 



410 

35 




5 24 

40 




197 

15 

200^/1 i ho •- _J 


644.5 cc/sec 

3 25 

5 28 

30 

25 


92012 , 

t 
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Then air is injected at various water i low— races. The 
frsquency so^ctruin (fraquency versus magnitude) at 3ach 
water £ lo'V’/ rate for various air-riot; rates are tabulated 
in Tables 2, 3, 4^5 and 6. The water i lew races are calcu- 
lated froui various pressures inside tne cnamber. 

The calculations of water flow rates from various 
pressures are done by DEC- 1090, Computer. The program 
and results are given. The problem of nubble formation 
and growth are iiscussed in Appendix-I. The calculations 
of radius of bubbles based on air-t low-rates are alsc^ 
cJiscussed in Appendix-I. 


TABLE 2 : Frequency Spectrum for 321.38 cc/sec. 

vifat^r Flow Raj^e 


Air Flow Rate 

With Bubbl 

as Frequency 

1 

1 Without 

Bubbles 


Hartz ; 

f 

,,, ^ - 1 - 

Magnitude 

'iFraq. ! 

Hertz ! 

Magnitude 


148 

50 mV 



128.9 cc/sec 

158 

187 

9 2.5 mV 

75 mV 

158 

100 mV 


261] to 

30 6j 

75 mV 

200 

80 mV 


363 75 mV 

392 9 2.5 mV 


257,8 cc/sec 397 200 mV 

552 50 mV 551 60 mV 


644' 

1 to 

30 

mV 

671 

f 



364' 

397. 

; to 

85 

mV 


386,7 cc/sec 552 50 mV 552 70 mV 

635) to 60 mV 

67^2 



367' 

394j 

! to 

100 mV 

515,6 cc/s 3 G 

650' 

670j 

1 to 

f 

40 mV 


548 50 mV 552 50 mV 




TABLE 3 


Fr>_qu3ncy Spectrum for 352.61 cc/sec 
water tlo’.r-rate 


Air Flow Rata 

J 

1 

1 u'ith Bubbl 

es Frequency 

Witinout 

Bubbl as 

Hertz 1 K 

1 1 

t f 

agnitude 

Freq. ; 

M irtz 1 

Magnitude 


391 to 

400 

30 

mV 

150 

40 mV 

128.9 cc/sec 

46^ to 

476 

25 

mV 




646 to 

756 

60 

mV 




360 to 

408 

65 

mV 

300 

50 mV 

257,8 cc/sec 

450 to 

475 

40 

mV 




648 to 

750 

60 

mV 




348 to 

454 

70 

mV 

300 

60 mV 

386,7 cc/sec 

640 


30 

mV 




630 


30 

mV 




764 


40 

mV 




343 to 

480 

70 

mV 

305 

55 mV 

515.6 cc/sec 

640 


45 

mV 




630 


40 

mV 




7 65 


50 

mV 




% 
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TABLE 4 ; Fr.^qu.Bncy Spectrum for 380. S6 cc/sec. of 

wa.t-3r flow-rate 


1 

I 

Air I low Rate 

r 

i 

1 

i 

I 

■ Vi tin 

Bubbl.as FroquonGy 

1 

1 

! without 

1 

Bubbles 

Hnrtz 

; Magnitude 
( 

( 

jFr-aj. ; 

iHertz 1 

Magnitud a 


195 

80 mV 




355 

40 mV 



128.9 Gc/s< 3 c 

405 

30 mV 

441 

30 mV 


512 

80 mV 




590 

SO mV 




350 

30 mV 




150 

30 mV 




197 

20 mV 



257.8 Gc/soc 

40 2 

45 mV 




470 

30 mV 

436 

35 mV 


575 

30 mV 




852 

50 mV 





136 

60 mV 




19 2 

40 mV 




347 

80 mV 



386.7 cg/ssg 

507 

35 mV 

438 

35 mV 


567 

35 mV 




850 

40 mV 




87 2 

40 mV 




142 

45 mV 



197 

30 mV 



332 

85 mV 




30 mV 

‘nc htO 

45 mV 

515.6 GG/saG 5 qq 

30 mV 



550 to 57 5 

35 mV 



847 

40 mV 



865 

40 mV 
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TABLE 5 : Frequency Spectrum for 407.15 cc/sec of 

water f low-rate 


Air F lo'!'7 Rate 

i 

i 

f 

j ;ith Bubbles 

1 

Frequaecy 

t 

I 

t 

J Without 

Bubbles 


1 Hertz J 

J f 

- . L 1 

Magnicude 

!Freq. 

iHertz 

Magnitude 

t 

i 


212 

25 mV 




385 to 437 

25 mV 

39 2 

30 mV 

128.9 cc/sec 

437 to 450 

600 to 765 

20 mV 

80 m.V 

552 

80 mV 


•335 to 900 

80 mV 




1097 to 1125 

40 m.V 




257,8 

cc/sec 

370 

650 

870 

1097 

216 

to 432 
to 760 
to 907 
to 1142 

20 mV 

30 mV 

65 mV 

75 mV 

30 mV 

393 

551 

30 

70 

mV 

mV 




210 

30 mV 






37 2 

to 440 

30 mV 

394 

30 

mV 

386.7 

cc/sec 

655 

to 7 55 

60 mV 

553 

65 

mV 



870 

to 908 

7 5 mV 






1097 

to 1125 

30 mV 





210 

25 mV 



370 to 430 

50 mV 

390 

30 mV 

515.6 cc/sec 645 to 760 

70 mV 

55 2 

70 mV 

870 to 910 

70 mV 



1085 to IBO 

30 mV 
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TABLE 6 t Frequanc:^ Spectrum for 431.85 cc/sec of 
water flow— rate 


Air Flow P.^ice 

With Bubbles Frequency 

1 

1 

1 

\ //ithouc Bubbles 

I 

Hertz 

Magnitude 

1 

! Freq. 

! Hertz 

1 

J Magnitude 


149 

80 mV 




25 5 

40 mV 




257 

40 mV 



128.9 cc/sec 

400 

30 mV 




530 

50 mV 

565 

30 mV 


642 

40 mV 




860 

30 mV 

880 

30 mV 


1020 

20 mV 




257.8 

cc/s ec 

147.5 

240 to 260 

390 to 400 

550 to 570 

620 to 640 

850 to 860 

1020 

70 mV 

30 mV 

50 mV 

40 mV 

30 mV 

40 mV 

30 mV 

558 

88 2 

35 

45 

mV 

mV 



152. 5 

60 mV 






347.5 

40 mV 






380 to 400 

50 mV 




386.7 

cc/sec 

550 to 370 

30 mV 

555 

50 

mV 



620 to 640 

30 mV 






840 to 860 

40 mV 

887 

35 

mV 



1020 

20 mV 






155 

50 mV 






350 

50 mV 






387.5 

50 mV 






437.5 

40 mV 




515.6 

cc/sec 

530.5 

40 mV 

550 

50 

mV 



562 

40 mV 






635 

40 mV 






840 to 860 

50 mV 






1014 

30 mV 

886 

45 

mV 






CHAPTER IX 


RBSJLTS AND DlSCUSSIOr^S 
* 

% 

In the previous chapter. Table 1 gives frequency 
spectrum for various air flow-rates in stationary water. 

It is found that the frequencies did not vary much with 
various air-flow rates, 

Strasberg [.4] nas snown that for air bubbles in 
water at atmosphere fr = 330 cm/sec (1) 

Venera f — frequincy, r - radius of bubble. 

So for various radii of bubbles, the frequencies should 
vary as in eqn.C 1) . 

But from Table 1 for various air-flow- rates, the 
frequencies did not vary. This indicated that there is no 
change in radius of bubbles. 

For various water-flow— rates (from Tables 2, 3,4, 5 and 6) 
also, the frequ-Bncies did not vary with air-flow-rates. 

Furthermore data of Tables 2, 3, 4, 5 and 6 reveals that 
the frequencies are shifted with water flow-rates, 

Th a shifted frequencies for 3ach water flow- rates 
witn and without bubbles are given in Table 7, 
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F roiTi the water flow— rates, the velocity of water in 
the chamber can be calculated. The computer program with 
results are attached, 

A graph of frequency versus 4th power of velocity of 
water is drawn (Fig. 3) . It is found that in both cases 
they are lin" ar. 

It shows that the frequencies are shifting with 4th power 
of velocity of water. 

That means the shifted frequencies are proportional 
to the 4th oower of velocity of water. 

So frequencies due to flow, i.e., without bubbles and 
frequencies dua to bubbles shift to the 4th power of velocity 
of water 

I.e. f 0^ v 

The curve fitting calculations between flow-rate 

A 

and ( velocity; for the datas in Table 7 are computed. The 
computer program with results are attached. 



CHAPTER X 


COHCLUSIOISIS i^LJD RECOMidEHDATIONS 


10.1 Conclus Ions 

From this experiment it. is observed tnat 

( 1) The noise cechniques can be used to obtain 
information about nucleate boiling. It is snown that 
the noise techniques can be used to determine the frequency 
response (spectral density) of bubbles in nucleate boiling 
for various ± low-races of coolant. 

From the experimienca 1 observations it is found that 

(1) The frequencies due to bubbles/ frequency due to 
flov>/ shirts to the 4th power of velocity of liquid. 

(2) Some of th s p ?ak;s at high ?r flow rates are due to 
frequency shift of the peaks present at lovj flow rates, 

10. 2 Recommend ations 

In this experiment frequencies due to bubble collapse 
are found. When the liquid is flowing all these frequencies 
are shifting to its 4th power of Vc:locity of liquid. The 
dependence of shifting of frequency as a function of surface 
tension, viscosity and density is done by dimensional 
analysis which is discussed in Appendix II. The theoretical 
explanations '.are not yet studied. In future it can be 


studied 
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appendix I 


Radius of Bubble Calculation 

The problem of bubole formation and growth are 
discussed here to calculate the radius of bubble [14]. 

Figure (4. a) shows the bubble formation in a cavity. 
Let R ba the radius of cavity. Let 'r' to be the radius 
of bubble. 

0 be the angle of contact between water and glass. 

S be the surface tension of water. 

oe tne pressure of the liquid - dynas/cm 

2 

be the pressure inside the bubble - dynes/cm 

In Fig.(4.b) the bubble is free floating X'/ithin a 
continuous liquid phase. The forces acting on that bubble 
are 

( i) The force dua to the pressur3 of the gas 

2 2 
(dynes/cin ) inside the bubble = nr P .. 

g 

( ii) Th 2 force due to the pressure of the liquid 

2 

(dynes/cm ) acting on the bubble = n P^ . 

(iil) The surface tension force S (dynes/cm) on the 
liquid-gas interfaces of the bubble = 27 cRS cos0 • 

Considering one half of tne bubbl a balance of 
the forces acting on it is represented by 
2 2 

■pcrP = nr + 2%BS cosQ 

2 27i:RS cos© 

r = “t; — — s — 


( 1 ) 


FIGURE 4(a) t BUBBLE FORMATION 



S coaO 


FIGURE 4(b) t FORCES Ca< A FREE BUBBLE WIIH IN A LIQUID 
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In 3qn.(l)/ ?., 3, Q are known. 
calculated rrom a: r-r lov7-rate. Therstore radius of the 
oubbl'=‘ can ca leu la tad. 

I'JOW/ 


Sjccess pressure 
inside the bubble 
i.e. (P -Pf) 


. f lo^v-rate 2 

^ area of cross-section 
of air rlow 


Jensity or air 


Substituting (2) in (1) radius of the bubble is calculated. 
The calculations of radius of bubbles based on air-f low- 
rates are done in D3C-1090 computer. Tne computer program 
with results are attached. 



APPTSi'.'OIX II 


It is tirst tnought ttie saitting of frequencies may 
oe also as a tunction of surface tension, density, in 
addition witn velocity of water. 7in atteir.pt is made to 
explain tnis based on dirrensional analysis (12) and (13), 

Frequency = f (P, s, v) 

f = k v^ ( i: 

v-velocity has dimension (LT”^) 

— 3 

P-density is mass per unit volume, (xiL ) 

S-surface tension (i.e. Force Per Unit Length) 

- 2 —2 

wnich has dimension MLT /L = MT 


(MT”^)^ (LT“^)^ 


-1 = —2b — c 

0 = a + b 

O = —3a c 


( 2 ) 

(3) 

(4) 


Solving eqns.(2)/ (3) and (4) it is obtained that a = 1, b = -1, 
c = +3. 

-1 3 

Sc eqn.(l) becomes f = kpS v 

^ , 3 

f = k p v 


By introducing one more parameter viscosity, in addition 

with surface, tension, density and velocity of the water, by 

^ n 4 

jL/ 

dimensional analysis it obeys the formula f = s— = 

s 

where p viscosity of the liquid. 



* PRUGRa.H T.l calculate the ELOWPAtE AMD VELHCITY of i^ATER From the * 

* MANOMETER HEIGHT f 

*♦*♦♦*♦ ?*♦*♦♦♦*!*♦ »#♦»»»♦♦♦«:♦»♦*♦»♦♦♦*♦»♦♦♦♦♦♦♦♦♦ 

* TO find the area op CROSS-SErTTDf! nF THE PIPE * 

PADitJS = 0.63 

Pi = 4,0 * A7ANC1.0) 

area = PT*RADITJS*RADIUS 

* TO calculte the flow-rate and velocity of water * 

*^:^tilc*t*t*****l¥***t*^***^**^^***^*ili*:^*t>¥******m:**3l::$:^)tf)ilft*t*1li*****'^******* 

0 = 980,0 
DHG = 13,6 
wRTrE(35,20) 

20 FbRMATC//////////////,l4 X, 'CALCUI.AtION OF FLOW-RATE AND VELOCITY 

1 OF waters/) 

WRirE(35,30) 

30 FOPMATf 10X,60C *♦')) 

MRirE(35,25) 

25 FORMATdOX, '♦ PRESSURE * FLOW-PATE » VEL 

lO^lTY *S 
write (35r30) 

DO 10 T = 5/9 
height = FLOAT (I) 

FRATE bSORTCHEIGHT *area*area*g»dhg) 

VH20 = FRATE/78,54 

WRITE (35»50)HEIGHT,FRATE/VH20 

50 F0RMAT(10X,'*'.F8,2/2X,»CM'.0F Hg. S 3X , SF8 ,2 , ' CC/SECS 

i 3Xf '*SF7,2, " CH/SEC *') 

10 CONTINUE 

WRTTEC35,30) 

STOP 


end 


w f ‘t 1 o / n ' -f * j 


i ii4567B9012 34 5 W« 01JTTJ^5 “6''7 3 Q U 1 2 4' S 'b 


CAr.iCnLaTTOM OF flow-rate AND VELOCITY OF WATER 



PRESSURE 


% 

FLOW-RATE 


velocity 





5'. 00 

CM.nF 

«a. 

% 

321.88 

cr/SEc 


4.10 CM/SEC 



6.00 

C*<.OF 

Hqr, 

% 

352.61 

CC/SEC 


4.49 CM/SEC 


♦ 

7'. 00 

CM. or 

Hg. 


380.86 

CC/SEC 

Jff 

4.85 CM/SEC 



8.00 

c^.or 

Hq. 

$ 

407.15 

CC/SEC 


5.18 Cm/SEC 

♦ 


9.00 

CM. OF 

Hq, 

♦ 

431 .85 

CC/SEC 

♦ 

5.50 CM/SEC 





iOil 

idi'j' 

i03i'-' 

(050 ' 

leso* 

• '' 

fl3Q' 

»140" 

»16U' 

1171^' 

H8W*' 

lioi#*' 

J20&'' 

j220‘' 
>23*^“ 
12 4 &' 
|25a‘r' 
| 260 "' 
1270 ?" 


**i^t****************t*t*****'^***'$**^ik*******y^^i^t***t***¥********t**f**** 
* PROGRAM TO FIND THE RADIUS OF BUPBI.ES FOR VARIOUS AIR FLOWRATES * 

**iif*****t ************* t**************^**t^H!^t ********* t*t****tt***t***** 
!4RTTE(22,25) 

25 PORHATC2X# // // // // // // // / ) 

WRTrE(22,10) 

10 FOP'^aTf9X# 'CALCULATED RADIUS UF BURBLES FOP AIR FLOwPATES ' , /3 

¥RTrE(22,30) 

30 POPMATC8X»47C '♦')) 

■4RirE(22,3S) 

35 FDRMATCSX,'^ air FLOWRATE * PADlffS OF BURBLES *' 5 

wRTXE( 22,30) 
no 20 T s 1,5 
P a PLOATCD* 128,9 
p =0.25 
T = 75.0 

COSTRE = C0SD(8.5) 
o2 = 2.0»RfT*CoSTHE/(0.0254*F»F) 
radius = S0B7(p2) 
fiJRTTEC22,50)F, RADIUS 

50 POPMATCSX, '♦'F!0.2,2X, 'CC/SET »',F15,5,' CM. *' ) 

20 CONTINUE 

MRTTEC22.30) 

STOP 

FND 


I' 


I 


C/iLrUTJiTED RADIOS OF BUBBLES FOP RiP FLOWRATES 


* ATR Flowrate ♦ radius of bubbles * 


♦ 

128.90 

cc/sec 

f 

0.2964S 

CM. 



2S7’.8n 

CC/SGC 

3jC 

0. j4822 

cm’. 

* 


386'.70 

CC/SEC 


0,09R82 

CM. 


t 

515'.60 

CC/SEC 


0, 07411 

CM. 


* 

b44’.50 

CC/SEC 

#3 

0,05029 

CM. 

3|t 


*'¥^*******^ ************** t******^***** 




\2i> 
>3C 

>5 J 
>6a 

n%>' 

^ 8 w^ 

f9w' 

1 0 ‘*» * 

• %'**' 

■ 2^' 
3 'j 
4w' 
5w' 
6 w 

7irf 

8W*' 

9t#t’ 

Owf" 

Jt#-> 

2wf: 

3v‘ 

4u’ 

:5W»~ 

6w 

7^"’ 

8wt" 

9«#‘ 

lOyr 

120'" 

(3W^ 
t40f' 
150 '.t 

t7W^ 
1 80.'^ 
I90('I» 

lO«€» 

\lQt 

^2Cli 

i30i 

I4©f4 


* THIS PROGRAM Is DESIGNED Tn FIND THE LEAST SQUARES PDLYNnMTAL * 

* THAI adequately REPRESENT A GIVEN SET OF DATA, THIS »itLL FIT ♦ 

* LEAST SQUARES POLYNOMIAL UPTn AND INCLUDING DEGREE lO * 

****************** *****t**V********* ******* ff***$*^**********t****** 
DIMENSION A(11,12),XPX(U,11),XPYC113 ,XC20)»Y(20),BC11) 
hpEN(UNIT=3l, ACCESS = 'SEOlN ' , FILE= *EOR23 . DAT ' ) 

DPENCUMIT = 5, DEVICES 'DSK',FILE = 'F0R33.DAT' ,ACCESS='SEQnUT' ) 
DO 121 INT = 1,9 
1 READC31 ,*)NUM 

t^-itl*^***^:*************************^* ********************************** 

* HEAD the OBSERVATIONS AND CALCULATE THE COEFFICIENTS FOR THE * 

* uDp-^AL EQUATIONS TO FIT A FIRST OFGPEE POLYNOMIAL * 

********************************************************************** 

SUMY = 0.0 
SUMX = 0,0 
suMXy = 0.0 
SUMXS = 0.0 

READOI,*) C XCl),y(I),I=l,NUM3 

DO 4 I =1,MUM 

type ♦, X(I),Y(I3 

SUMY = SUMY + YCI3 

SUMX = SUMX + X(I) 

SUMXY = SUMXY + XCT)*YCI5 
SUMXS = SUMXS + X(I3*X(I) 

4 RMSSQ = 10.0 E 08 

K =l 

N = K + 1 
M = M T 1 
XPX(1#1) = NUM 
XPXCl,2) = SUMX 
XPX(2,i) = SUMX 
XPX(2#2) = SUMXS 
XPYCl) = SUMY 
XPY(2) = SUMXY 

******************************************************************** 

* transfer IHE XPX and XPY QUANTITIES TO THE MATRIX AND SAVE THE ♦ 

* XPX and XPy QUANTITIES. THE GAUSs-JORDAN METHOD WILL PERFORM » 

* elementary transformations until the solution of the normal * 

* EQUATiblNS IS IN A(I,M) , 1=1, N, * 

******************************************************************** 


'; 49 v ' 

:i50m 

:<5 1 V' 

^52‘*" 

j53u‘ 

i54«jc 

s 55 w 

Si 56 w 

?57y'' 

■58v' 

!59v' ■ 

it60af 
|6 1 w.' 
(62W’ 
|63U* 
i64C»;" 
1655#'^ 

i660‘" 

i 67 ^iT 

i6 8 w*‘ 

169^' 

|70C-' 

!7iy' 

>74y*‘ 

»75C#'" 

1760'" 

>7707. 

|78y7 

>7907 

|8O07 

iSiy 7 

|82fe7 

|83§7 

p 40 # 

18507 

■8607 

18707 

|88@7 

■BOCli.. 


IB DO 12 I = l^N 

DC ID J s 1 ,M 
to aci.J) = XPX(l,J) 

12 = XPYCI) 

^4***^*****^'¥^***^***^*****t!l!******)¥*^ilf**********it^**-******t*******Jlf^ 

* subroutine calling to solve the normal equations to Get the » 

* Least squares estimates of the parameters =f 

CALL GSGnK(N,M,A) 

SUMSO = 0.0 
no 16 I 3 1 ,NUM 

prod = 0.0 

no 17 J 3 1 ,K 

17 PROD = PR 0 D+A(J+ 1 /M)*XCI)»fa 

YHAT = A( 1 ,M)+PR 0 D 

16 SUMSO = SUMSQ+(Y(I)-YHAT )**2 

TF(NUM-K- 1 .E 0.03 go to 42 

*t*******¥**************t*^*^******t********************^**************** 

* calculate the current modified sum of squares * 

»♦♦♦*♦♦♦♦ 

CMSSO = SUHSQ/{NUM-K-l) 

* IF THE CURRENT MODIFIED SUM OF SQUARES IS GREATER THAN THE PREVIOUS * 

* modified Sum of squares, an adequate degree polynomial has been ♦ 

* obtained, so go to 42 ; Otherwise, WRITE the degree of the polynomial * 

* .TUsT Obtained, the corresponding least squares estimates of the * 

* parameters, AND THE CURRENT MODIFIED SUM OF SQUARES ♦ 

IF(CMSSQ.GE.PMSSO) GO TO 42 
WRITE C 5 , 60 )K 

60 F 0 RMAT(/ 2 X,'THE COEFFICIENT OF THE LEAST SQUARES POLyNiMTAL 

I OF DEGREE%I 2 ,' are',/) 

00 15 I 3 1 ,H 
IM = I-l 

15 WRirEC 5 , 65 )IM,ACl,M) 

65 F 0 RMATC 2 X,'BETA(',I 2 ,')=',F 12 . 6 ) 

WRirE( 5 ,? 0 } CMSSO 

70 FORMAT!/, 2 X, 'CMSSQ= ' ,F 15 , 4 ,//) 

PMSSQ 3 CMSSO 
IFCK.EQ.IO) GO TO 1 


^ tt ****** ************************************************ ****’¥ t *1 f ******^ 

* INCREASE IHE DEGREE OF POLYNOMIAT, BEINC FITTED ,K ,B!rlt,D IRE NE»tf * 

* put.iromial E0UATI0NS,AMD then return id statement 18 ♦ 

^^ilf:^ilf.ltf*$llf*^f*t*************t*******!t:************************************* 

K = K + 1 

N = K +1 

M = N + 1 

Ypi = K + 1 

fCHl = fC - 1 

DO 20 1= l,KMl 

KPXCKPl,!) = XPX(K,I+1) 

20 V?X(I,KP1) = XPX(KP1,I) 

SUMl = 0.0 

SUM2 s 0,0 
SUM3 = 0.0 
DO 21 I = 1,NUM 
XK = X(I)»*K 

SUMl = SUMl + XK*X(I)*»CK-n 
SUM 2 = SUM2+XK*XK 

21 SUM3 = SUM3 + XK»XCI) 

XPXCKPl.K) = SUMl 
XPX(KP1,KP1) = SUf42 
XPXCK,XP1) = SUMl 
XPyCKPl) = SUM3 

nil 27J = 1,K 
27 R(J)= A(J,M) 

CO TO 18 

If ^^ itiiltiif ^ llt *^;^^***********************^^************************************** 

* PRINT VEl<aCITY**4,FRE0UENCY(EXPERTMENTALJ ,CALCULATEDfERROR FOR THE * 

* POUINOMIAO ADEOUATEUY REPRESENTING DATA * 

************************ ***1F************************************** 

42 XMt = K-1 

NRirE(5,80)KMl 

SO F0RMAT(3X,//,3X» 'THE POLYNOMTAt OF DEGREE ' ,13 a X, ' ADEQUATELY 

t represents data') 

HRITE(5,85) 

85 FbRMATC//4X, '(CC/SEC)**4'6 x» 'FREQUENCY', 6X, 'CALCULATED 3X, 

13X» 'ERROR'/) 

DO 25 I = 1,NUM 
prod « 0.0 
no 26 Js l,KMl 




i36y 



137-' 




26 

PROD = PROD + B(J+1 )*X(I)*»J 

t39u' 


Yri^r = S(l) + PROD 

i40U* 


DIFF = 5fCI) - YHAT 


25 

■AmrrE(5,90)XCI),Y(I) ,yHAT,DlFF 


QO 

fORHATC4F15,6) 

t43ir 

I2t 

eonitnue 

!44w" 


SfnP 

f 


fnd 

i46w' 


subroutine GSG0R(N,M,A) 

!47v' 


niME»rsiONAciirr2) 

Hsitr 


DO 16 K = 1,N 

i49fr 


ft 

|50U' 


DO 9 a = KPl,M 



A(K,.I) = ACK,J)/ACK,K3 

t52U^* 

9 

CONTINUE 

jSBm ■ 


DO 12 T = 1»N 



TF ( I .EQ.K) GO TO 12 



no 14 J = KP1,M 

t56&f 


&{I,J) = A(I,J)- A(I.K)*A(K,J) 


14 

CONTINUE 

i58u' 

12 

CONTINUE 

i59w' 

1 b 

continue 

ibOtJf' 

1 i 

RETURN 



END 


THE CnE«-FTCIEHr OF THE LEAST SQUARES PDLYnIMIAL OF DEGREE 1 ARE 

tiETtC 0)= 452,147750 

BETAC 1)= 0.737672 


CmSSQs 


8g.6503 


THE POLYhOmIAL of DEGREE i ADEQUATELY REPRESENTS DATA 


(CC/SEC7»*4 

282,540000 

406.430000 

553.310000 

915.060000 


rREQUE^fCY 

660.000000 

760.000000 

850.000000 

mo.oooooo 


calculated 

660.569560 

751.959720 

860.308960 

1127.161800 


ERROR 

-0.569557 

8.040276 

•10,308960 

2.S38242 


THE COEPFICTENT OF THE LEAST SQUARES POLYnIMIAL OF DEGREE 1 ARE 


bexac 0)= 
BETaC i)= 

CMSSQs 


319.751980 

0.772799 

1 8 , 4l 03 


TriE POlYMOMIAL of degree 1 adequately represents DATA 


(CC/SEC)**4 

282.540000 

406.430000 

719.980000 

915.060000 


FREQUENCY 

540.000000 

630.000000 

880.000000 

1025.000000 


calculated 

538.098540 

633.840580 

876.151640 

1026.909200 


ERROR 

1.901459 

-3,840584 

3,848358 

-1,909225 


The coefficient of the least squares PQLYnIMIAL of degree 1 ARE 


BETaC 01= 
BETAC 11= 

CMSSQs 


134.685480 

0.784262 

30.07J5 


The coefficient of the least squares POLYNIMIAL of degree 2 ARE 


BETaC 01= 113.239960 
BETAC i)= 0.867535 
BeTaC 2)= -0.000069 


CMSSQ: 


24.1107 


THE coefficient OF THE LEAST SQUARES PULYNiMlAL OF DEGREE 3 ARE 


BETAC 

BRTAC 

BRTaC 

BETAC 

CMSSQs 


01 = 
11 = 
21 = 
31 = 


52.615763 
1.232747 
-0.000732 
0. 000000 

19.2016 


the polynomial of DEGREE 3 ADEQUATELY REPRESENTS DATA 
CCC/SEC1**4 FREQUENCY CALCULATED ERROR 


282.540000 

406.430000 

553.310000 

719.980000 

915,060000 


460.000000 

570.000000 

700.000000 

850.000000 


350.790620 

457.47U10 

573.060750 

698.354230 

850.323300 


-0.790615 

2.528893 

-31060753 

1.645767 

-0,323296 


THE COEFFICIENT OF THE LEAST SQUARES POLYNlMlAL OF DEGREE 1 APE 


CMSsas 


2000.6265 


rR£ CnKFFTCTE’^r UF THE least squares POLYNImIAL of degree 2 APE 


betac 0)= 

BET«i( Ds 
BE'i'A C 2 ) s 

CMSSQs 


77.955549 
0. 123715 
0.000512 

1787.4162 


the pot.ymomial of degree 2 adequately represents data 


(CC/SEC)**4 

282.540000 

553.310000 

719.980000 

915.090000 


frequency 

148.000000 

330.000000 

402.000000 

630.000000 


calculated 

153.819940 

303.2919O0 

432.656260 

620.231900 


ERROR 

-5,819944 

26,708099 

-30.656258 

9,768097 


THE COEFFTCIENT OF THE LEAST SQUARES POLYNIMIAL OF DEGREE 1 APE 

BETAC 0)= -308.134980 
BETA( 1)= 0,919142 


CmSSQ= 


61.0085 


the polynomial of DEGREE 1 ADEQUATELY REPRESENTS DATA 


(CC/SEC)»»4 

553.310000 

719.980000 

915.060000 


frequency 

197.000000 

360.000000 

530.000000 


calculated 

200.435610 

353.629060 

532.935330 


ERROR 

-3,435612 

6,370945 

-2.935326 


the pulynomial of degree 0 adequately represents data 


CCC/SEC)*#4 FREQUENCY CALCULATED ERROR 

719.980000 200.000000 353.629060 -153.629060 

915.060000 380.000000 532.935330 -152,935330 

The COEFFtCTEMT of the least squares POLYnImIAL of degree 1 APE 

BEIAC 01= -29,609509 

BETAC 11= 0.817648 

CmSsQs 237.4820 


THE polynomial OF DEGREE 1 ADEQUATELY REPRESENTS DATA 


(CC/SEC)*f4 

282,540000 

406,430000 

553.310000 

719.980000 


frequency 
goo, 000000 

295.000000 

441 .000000 

550.000000 


calculated 

201.408770 

302.707180 

422.803320 

559.080730 


ERROR 

-1.408768 

-7,707180 

18.196678 

-9.080727 


THE COEFFTCTENT OF THE LEAST SQUARES POLYNIMIAL OF DEGREE 1 ARE 


betac 0)= 

6ETAC l) = 
CMSSQs 


172.022330 

0.787152 

343,3217 


THE coefficient OF THE LEAST SQUARES POLYNIMIAL OF DEGREE 2 ARE 
Betac oi= 8i, 285402 

r 1 1= < < HQAaA 


TriE PULYMOMlAii OF DEGREE 2 ADEyUATELY REPRESENTS DATA 


CCC/SFC)=‘‘*4 


FREQtJEMCY 


CAtiCULATED 


ERROR 


782.540000 

406.430000 

553.310000 

719.080000 

915.060000 


3R0. 000000 

500.000000 

610.000000 
76U. 000000 

875.000000 


379.787900 

495.887630 

621.850260 

749.434910 

878.039280 


0.212101 

4.112370 

-11,850758 

10.565086 

-3.039276 


THE CnEFfiCTENT OF THE LEAST SQUARES POLYNlMlAL OF DEGREE 1 


BETaC 0)= 
BETAC 1)= 

C^SSQ=: 


•170.478360 

0,784452 

28.3857 


TriF polynomial of degree 1 adequately represehts data 


CCC/SEC)4*4 

406.430000 

719.980000 

915.060000 


FREQUENCY 

150.000000 

390.000000 

550.000000 


calculated 

148.346440 

394.311340 

547.342220 


ERROR 

1,653564 

- 4:311337 

2.657776 


ARE 



KieTp-nSt-M- <Rl-I>£T 



